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ABSTRACT 
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We report on a 40 ks long, uninterrupted X-ray observation of the candidate supergiant 
fast X-ray transient (SFXTs) IGR 116418^532 performed with XMM-Newton on Febru- 
ary 23, 2011. This high mass X-ray binary lies in the direction of the Norma arm, at an 
estimated distance of 13 kpc. During the observation, the source showed strong variability 
exceeding two orders of magnitudes, never observed before from this source. Its X-ray flux 
varied in the range from ^0.1 counts s^^ to ^15 counts s^^, with several bright flares of 
different durations (from a few hundreds to a few thousands seconds) and sometimes with 
a quasi-periodic behaviour This finding supports the previous suggestion that IGR J16418- 
4532 is a member of the SFXTs class. In our new observation we measured a pulse period of 
1212±6 s, thus confirming that this binary contains a slowly rotating neutron star During the 
periods of low luminosity the source spectrum is softer and more absorbed than during the 
flares. A soft excess is present below 2 keV in the cumulative flares spectrum, possibly due 
to ionized wind material at a distance similar to the neutron star accretion radius. The kind of 
X-ray variability displayed by IGR 116418^532, its dynamic range and time scale, together 
with the sporadic presence of quasi-periodic flaring, all are suggestive of a transitional accre- 
tion regime between pure wind accretion and full Roche lobe overflow. We discuss here for 
the first time this hypothesis to explain the behaviour of IGR J16418^532 and, possibly, of 
other SFXTs with short orbital periods. 
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1 INTRODUCTION 

IGR J16418-4532 is a transient X-ray source discovered in the di- 
rection of the Norma Galactic arm during INTEGRAL obser - 
vations performed on 2003 February 1-5 iXomsick et al J l2004h . 
Its propo sed optical/infrared counterpart (2M ASS J16415078- 
4532253. Ictiatv et al.|[2003, iRahoui et al.ll2008l) . positionally con- 
sistent with the accurate localization obtained with the Swift X- 
Ray Telescope ^Romano et al.ll2oTlh . is indicative of a high mass 
X-ray binary (HMXB) at a distance of about 13 kpc. The short time 
interval of X-ray activity observed with INTEGRAL suggested 
to class ify IGRJ16418-45 32 as a supergiant fast X~ray transient 
(SF XT,ISguera et al.ll2006h . 

IWalter et al.l ( l2006h reported the presence of pulsations at 
1246±100 s in an XMM-Newton observation obtained in 2004, 
while a periodicity of 3.753±0.0q4 days ICor bet et al. 2006) or 
3.7389±0.0004 days dCorbet et alj|2006l. iLevine et al.ll201 10 was 
obtained from an analysis of the long term light curves from the 
Swift Burst Alert Telescope and the RXTE All Sky Monitor 



data (ASM), respectively. This modulation has been interpreted as 
the orbital perioqj of the binary system. 

Following the MAXI satellite discovery of renewed X-ray ac- 
tivity from the region of IGR J16418-4532 in February 2011 (as 
reported in lRomano et alj|201 10 . we performed a Target of Oppor- 
tunity XMM-Newton observation. The new XMM-Newton 
data show variability spanning an unprecedented range in this 
source, thus supporting its SFXT nature. Furthermore, the frequent 
presence of strong quasi-periodic flares suggests that this long spin- 
period pulsar in a relatively compact HMXB (as implied by the 
short orbital period), might be in an intermediate accretion regime 
between Roche-lobe overflow and wind accretion. 



2 OBSERVATIONS AND DATA REDUCTION 

IGRJ16418-4532 was observed with XMM-Newton between 
2011 February 23 (at 13:55 UT) and February 24 (at 00:32 
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the formal discrepancy between the two periods is probably due to an 
underestimation of the respecive errors. 



2 L. Sidoli, et al. 



UT), with a net exposure of about 39 ks. The X M M-Newton 
Observatory carries three 1500 cm^ X-ray telescopes, each 
with an European Photon Imaging Camera (EPIC) at the focus. 
Two of the EPIC u se Metal Oxide Semi-conductor (MPS) CCD s 
jTumer et ani200lh and one uses a pn CCD ( Struder et alllaOOll). 
Refle ction Grating Spectrometer (RGS) arrays tden Herder et alj 
l200lh are located behind two of the telescopes. 

Data were reprocessed using version 1 1 .0 of the Science Anal- 
ysis Software (SAS). Both MOS and pn operated in Full Frame 
mode and used the medium thickness filter. Extraction radii of 40" 
and 1' were used for the source events, respectively for the pn 
and MOS cameras. Background counts were obtained from simi- 
lar sized regions, offset from the source position, and in the same 
temporal intervals, when dealing with time selected analysis. The 
background (selected with PATTERN=0 and in the energy range 
10-12 keV in the pn) showed evidence of flaring activity only dur- 
ing the first 12 ks of the observation. Therefore, we removed the 
corresponding time interval in most of our analysis, as discussed in 
detail in the next section. 

Response and ancillary matrix files were generated using the 
SAS tasks RMFGEN and ARFGEN. Using the SAS task EPATPLOT, 
we found that MOS spectra were affected by pile-up. Thus, we re- 
port here only on EPIC pn spectroscopy (adding MOS data did not 
improve the spectral fitting), while for the timing analysis both the 
two MOS and the pn were considered. Spectra were selected using 
patterns from to 4 with the pn. 

To ensure applicability of the statistics, the net spectra were 
rebinned such that at least 30 counts per bin were present and such 
that the energy resolution was not over-sampled by more than a fac- 
tor 3. All spectral uncertainties and upper-limits are given at 90% 
confidence for one interesting parameter. In the spectral fitting we 
used the photoelectric abs orption model PHA BS in XSPEC with the 
interstellar abundances of I Wilms et alj ([2000l). The R GSs were op- 
erated in spectroscopy mode (Id en Herder et al .120011) , but given the 
high absorbing column density, they did not detect the source. 

To better investigate the long term behaviour of IGR J16418~ 
4532 we reanalysed also the 2004 XMM -NewtonfEPIC observa- 
tion, first reported bv lWalter et alj ( l2006l) , with the same procedures 
and selection criteria used for the 201 1 datfl 



3 ANALYSIS AND RESULTS 
3.1 Spectroscopy 

The EPIC pn, background subtracted, light curves of IGR J 164 18- 
4532 observed in 2004 and in 2011 are shown in Fig. [T] Letters 
mark time intervals showing different kinds of X-ray activity within 
each observation. In particular, strong flares are seen in time inter- 
vals B, D, and F, while in the remaining time intervals the source 
was in a low intensity state, with less frequent and fainter flares. As 
mentioned above, the first part of the 201 1 observation (time inter- 
val D) was affected by high background level, so in the following 
we will not use this time interval, except when extracting the spec- 
tra from the peaks of the bright flares (for which the background 
contribution is negligible). 

We first analysed the average spectrum of each observation, 
corresponding to a net integration time of 19.6 ks for the 2004 data, 

except for the pn source extraction region, for which we adopted a radius 
of 20" to avoid stray light contamination from a bright transient source 
located outside the field of view. 



Table 1. Spectral results of the time averaged spectra of the 2004 and 2011 
XM M-Newton observations. An absorbed power law model was used, 
r is the power law photon index. Flux is in the 1-10 keV energy range in 
" ^ and is corrected for the absorption, Njj (in 



units of 10^2 , 



")■ 



Parameter 


2004 (A+B+C) 


2011 (E+F) 


Nh 


15 ± 1 


7.4 ±0.3 


r 


1.3 ±0.1 


1.28 ±0.05 


Unabs. Flux 


1.8 


1.5 




1.147/166 


1.252/210 
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Figure 2. Spectral results of the average spectra extracted from the EPIC pn 
observations in 2004 {upper panel) and in 201 1 (lower panel). Count spec- 
tra are shown, together with the residuals (in units of standard deviations) 
of the data to the absorbed power law model (Table[T). 



and 23.8 ks for the 2011 data (intervals E and F). The best fit pa- 
rameters obtained with an absorbed power law model are listed in 
Table [T] and the spectra shown in Fig. [2] A similar spectral slope 
is seen in the two observations, while the absorption in 2011 was 
about a half of that observed in 2004. While the 2004 spectrum 
is reasonably well fitted by the power law model (although with 
marginal evidence for a soft excess below 2 keV), structured wave- 
like residuals over the entire energy range are evident in the 2011 
spectrum, with a positive excess below 2 keV more pronounced 
than in 2004. 
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Figure 1. EPIC pn, background-subtracted, light curves of IGRJ16418^532 in the 0.3-12 keV energy range. The bin size is 50 s and the time axis is in UTC 
hours of the dates indicated in the panels. Letters indicate the different intervals used for the time selected spectroscopy. Note that also during the non-flaring 
time intervals, the source is significantly detected, with average count rates 0.185±0.006 counts s~^ and O.331±0.005 counts s~^ in 2004 (interval C) and 
2011 (interval E), respectively (0.3-12 keV). 



We next considered temporal selected spectra, extracting dif- 
ferent spectra from the flaring intervals and from the low intensity 
emission, within each observation. 

We first considered the 2004 observation, when IGR J16418- 
4532 showed a bright flare with a complex shape (interval B, Fig.[T] 
left panel) in the middle of two low-intensity time intervals (A and 
C). We extracted a spectrum from each time interval (A, B and C), 
with the one of interval B formed by selecting only the events cor- 
responding to time intervals with count rate above 0.5 counts s^^. 
The results for these spectra, fitted with an absorbed power law, 
are reported in Table |2] We also extracted a further low intensity 
spectrum from time intervals A and C, but now excluding the two 
faint flares that occurred in time interval A. This "cleaned A-l-C" 
spectrum is indeed more representative of the low intensity emis- 
sion in 2004. A comparison with the flare spectrum (B) indicates 
that the IGR J16418-4532 X-ray emission is less absorbed during 
flares. There is also a possible hint for a flatter spectral slope when 
the source is brighter, but the uncertainties are too large to draw a 
firm conclusion, based on the 2004 observation only. 

We next considered the observation performed in 201 1. In or- 
der to investigate the spectral properties of the single flares we se- 
lected all time intervals with a count rate above I counts s~^ in the 
pn light curve (0.3-12 keV) binned at 256 s. This time selection 
yielded 16 single flares (6 flares in time interval D and 10 in F). 
An absorbed power-law gave almost always a good fit to each sin- 
gle flare spectrum. The best fit parameters of each flare are shown 
in Fig. |3] which shows no evidence for significant spectral vari- 
ability between the flares, even comparing flares in the first part 
of the observation (interval D) with those located in the second 
part (interval F). Most of the flare spectra require an absorption 
well in excess of t he total Galactic value towa rds IGR 116418^532 
(1.9x 10^^ cm~ MDickev & Lockmar]|l99(]|) . A black body model 
always results in a worse (or in an equally good) fit, giving densi- 
ties, on average, a half of the value obtained with a power law. The 
black body temperatures are ~2 keV, and the resulting black body 
radii have values of a few hundred meters at 13 kpc. 

We next extracted a cumulative "spectrum during flares" by 
summing all the spectra of the flares occurring in the time interval 
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Figure 3. Spectral results of the intensity selected spectroscopy of all flares 
observed in 2011. The assumed model is an absorbed power law. Upper 
panel shows the absorbing column density evolution with the flux corrected 
for the absorption, while the Lower panel displays the photon index versus 
unabsorbed flux (1-10 keV). Light triangles mark the spectral parameters 
of the flares extracted from the first part of the 201 1 observation (interval D 
in Fig. 1), while solid black squares indicate flares occuiTed in the second 
flaring part of the same observation (interval F in Fig. 1). 
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Table 2. Results of the time selected spectroscopy (letters mark the same time intervals displayed in Fig. 1). An absorbed power law model was used. Flux is 
in the 1-10 keV energy range in units of 10~^^ erg cm"^ s~^ and is corrected for the absorption, Njj (in units of lO^-^ cm"^). 



Parameter 


A 


B 


C 


A + C 
(cleaned) 


E 


F 


Nh 

r 

Unabs. Flux 


2Q+_l 

1 oc: + 0.25 
1.1 

1.035/41 


i.29«:i} 

6.5 

1.130/128 


0.7 

0.690/35 


19+2 

1 54+0-17 

^-^^^-0.22 

0.8 

1.085/60 


o 9+0.8 

0.65 
1.008/90 


-79+0.36 
°- '^-0.34 
1 -11+0.06 
-^■-^-^-0.06 
4.4 

1.196/176 



F, reaching a net exposure time of 5085 s. An absorbed power law 
gave a good fit above 2 keV, but showed some positive residuals 
(although never exceeding 3 o from the model) at softer energies 
(reduced x2=l-196 for 176 degrees of freedom, dof, see the upper 
panel in Fig. |4]l. These positive residuals suggest the presence of 
either (1) an additional soft component, (2) an ionized absorber or 
(3) a partial covering. Thus, we first added a black body compo- 
nent to the power law continuum, obtaining a better fit (reduced 
x2=0.987 for 174 dof; F-test probability of 5.53 x 10"*). This re- 
sulted in a softer and more absorbed power law, and a black body 
with temperature of 0. 17 keV and emitting radius of a few hundred 
km. 

The use of an ionized absorber (ABSORI in XSPEC: lDone et all 

in addition to a neutral one accounting for the interstellar 
absorption (PHABS), also resulted in a better fit with respect to the 
simple absorbed power law (reduced x^=0.950 for 174 dof; F-test 
probability of 1.99x 10^^). The absorber temperature was fixed to 
30 kK and the power law photon index of the photoionizing source 
was linked to the photon index of the continuum. The absorber ion- 
ization state ^ (defined as ^=hlvS? , where L is the X-ray luminos- 
ity of the ionizing source, n is the density of the absorber, and R is 
the distance of the ionized matter from the X-ray source), resulted 
in a value of 125 I45 erg cm s^^. 

A better fit with respect to the simple absorbed power law (re- 
duced x^=0.971 for 174 dof; F-test probability of 1.335x10"*) 
was obtained also using a power law continuum modified by 
a partial covering absorption (PCFABS in XSPEC), together with 
the usual absorption accounting for the interstellar value. Partial 
covering parameters are the equivalent hydrogen column density 
Nnpcfabs and the dimensionless covering fraction, f (0 < f < 1). 
All the best fit parameters of these models are summarized in Ta- 
ble [3] The counts spectra together with their residuals when fitting 
with three different continue are reported in Fig. |4] The residuals 
of the partial covering model are not shown because they look very 
similar to those obtained with the additional soft component (mid- 
dle panel in Fig.Q. There is no evidence of an iron emission line. 
Upper limits (95% confidence level) to the equivalent width (EW) 
of a narrow line from ionized iron (energy fixed at 6.7 keV) can be 
placed at EW<34 eV, while for a neutral iron emission line (energy 
fixed at 6.4 keV) at EW< 16 eV, assuming a power law continuum. 

The overall flares spectrum can be compared with the cumu- 
lative spectrum extracted from the 2011 low intensity state (time 
interval E, average count rate of 0.330±0.005 counts s"^), corre- 
sponding to a net EPIC pn exposure of 14.4 ks. An absorbed power 
law gave a good fit ixl = 1.008 for 90 dof) with the following 
parameters: ArH=(8.2lg;*) x 10^^ cm"^ r=1.57l[;;Jg. The aver- 
age observed flux was 3.7x10"^^ erg cm"^ s"^, while the flux 
corrected for the absorption was 6.5x10"^^ erg cm"^ s"^ 



Table 3. Spectral results for the cumulative spectrum during flares in time 
interval F (2011) fitted with models able to account for the soft excess. 
Flux is in the 1-10 keV energy range in units of lO^^^ erg cm^^ s^^ 
and is corrected only for the absorption Nh (PHABS in XSPEC). Absorbing 
column densities are in units of 10-^^ cm~^. 



Parameter 


POW + BB 


ABSORI * POW 


PCFABS * POW 


Nh 


9.41«-j« 


n 65+0-™ 

'-'•''0-0.50 


5 00+0.69 


^Habsori 




15.8+2-0 








125t«i 




Nnpcfabs 






Il.Tlt.? 


f 






0.59l°:l? 


r 


1 oq+0.08 
±.oo_g 


1 29+0-07 


1 42+0-12 
-'--^^-0.12 


kTbb (keV) 


o.i7lg:S^3 






Rbb (km) 


250l«o 






Unabs. Flux 


10.7 


3.1 


4.0 


X^/dof 


0.987/174 


0.950/174 


0.971/174 



(1-10 keV), translating into an average luminosity of 
1.2x10^^ erg s"^ at 13 kpc. Also an absorbed black body 
gave a good fit (xl = 0.924 for 90 dof), with iVn, of 
3.6±0.4xl022 cm~2, a temperature, kTbb, of 1.78±0.07 keV 
and a radius, Rbb, of 0.29±0.02 km (at 13 kpc). Finally, we fitted 
the low intensity spectrum with the best fit model obtained for 
the flares, i.e. an absorbed power law modified with a ionized 
absorber. This improved the fit quality (x^ ~ 0.794 for 88 dof, 
F-test probability of 2.75 x 10^'') and gave the following parame- 
ters: A'^H=(4±2)xl022 cmr'\ r=1.85 to fg, ionized absorption, 
-/Vnabsori, of (2.5 ^0 5) X lO^"* cm"^, ionization parameter, ^, of 
600 I335 erg cm s"^. If we fix the ionized absorption parameters 
A'nabsori and ^ values obtained for the cumulative flares spec- 
trum, we again obtain a good fit to the low intensity spectrum 
(X^ — 0.886 for 88 dof), resulting in a column density Nn=(2J 
+g g)xl0'^2 cm^2 and in a power law photon index r=1.7±0.1, 
consistent with previous results. Therefore, we conclude that the 
ionized absorber found in the overall flaring spectrum is consistent 
with being present also during the low intensity emission with 
similar properties, although the statistics does not permit to be very 
sensitive to its parameters. In any case, the low intensity emission 
is softer than the flaring spectrum. The low intensity spectrum 
fitted with these models is shown in Fig.|5] 

Finally, we also performed an intensity selected spectroscopy: 
a series of intensity selected spectra were produced, using intervals 
corresponding to EPIC pn count rates of <1, 1-5, 5-10, and >10 
counts s"^, when the data are accumulated with a binning of 256 s. 
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Figure 4. Cumulative X-ray spectrum of the flaring activity of the second 
part of the 2011 observation (flares in interval F). The upper panel shows 
the counts spectrum together with the residuals (in units of standard devia- 
tions) when fitting with an absorbed power law. The middle panel shows the 
residuals when fitting with a power law together with a black body. Lower 
panel displays the residuals when fitting with a power law continuum mod- 
ified by a ionized absorber (ABSORI in XSPEC, see Sect l3.1l for details). 

for both observations (2004 and 2011). The best-fit, obtained with 
an absorbed power law, is reported in Table |4l and it clearly in- 
dicates that the spectrum is harder when the source is brighter, as 
usually observed in accreting pulsars. 

3.2 Timing Analysis 

For the timing analysis we used both MOS and pn EPIC data. Ar- 
rival times were corrected to the Solar System barycenter, but not 
for the 3.7 days orbital motion, owing to the unknown parameters 
of the system. The observability of the periodicity at 1246±100 s 
dWalter et alj|2006l) is complicated by the presence of the strong 
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Figure 5. X-ray counts spectra and residuals (in units of standard devia- 
tions) for the low intensity emission observed in 2011 (interval E) fitted 
with an absorbed power law model (upper panel), an absorbed black body 
(middle panel), and an absorbed power law modified with a ionized ab- 
sorber (lower panel). 



flares which occur on a comparable time scale. Folding the data 
of the whole observation at different trial periods produces several 
peaks in the distribution. On the other hand, by restricting the 
period search to the central part of the 2011 observation (interval 
E), when the source was in a low intensity and non-flaring state, the 
X-ray pulsations are clearly detected at a period, P, of 1212±6 s 
(X^=252, for 9 dof). The corresponding folded light curve is plot- 
ted in Fig.|6] 

Also in the 2004 observation the source pulsations are more 
easily detected during the non flaring time interval. By a standard 
folding analysis of the time interval C we obtained a clear peak in 
the distribution at P = 1216±7 s (x^=257 for 9 dof). Based on 
this result, we could then refine the period estimate using the data 
from the whole observation, which gave P = 1213±2 s. The folded 
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Table 4. Results of the intensity selected spectroscopy. An absorbed power 
law model was used. Flux is in the 1-10 keV energy range in units of 



IGR J16418-4532 



10- 



erg cm 



of 10 cm ). Uncertainties on the unabsorbed fluxes are about 3%. 



Obs. 2004 


<1 c s^i 


1-5CS-1 


5-10 cs-i 


>10c s-i 


Parameter 











Parameter 



Nh 


o 9+0.8 
°-^-0.6 


O-»-0.3 


70+0.53 


3-O"-0.50 


r 


, ^7+0.13 
-0.10 


1 1 9+0.05 
-'-•-'-^-0.05 




n ss+0-09 


Unabs. Flux 


0.7 


5.0 


11.7 


20.0 




1.008/90 


1.147/209 


0.823/119 


1.190/137 



XMM 2004 



and is connected for the absorption, Njj (in units S ^ 









O 




Nh 


19+2 




"c 
E 




r 


1 t:4+0.17 
J^-^^-0.22 


1.29l«- 


o 




Unabs. Flux 


0.8 


6.5 




d 




1.085/60 


1.130/128 




in 



Folded period: 1216s 




Figure 7. IGRJ16418^532 pulse profiles at soft (0.3^ keV) and hard 
energies (4-12 keV) with XMM-Newton in 2004 (only the low inten- 
sity emission), together with their hardness ratio, obtained folding the light 
curve at a period of 1216 s. The zero phase is arbitrary. 
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Figure 6. IGR 116418^532 pulse profiles at soft (0.3^ keV) and hard 
energies (4-12 keV) with XMM-Newton in 2011 (only the low inten- 
sity emission), together with their hardness ratio, obtained folding the light 
curve at a period of 1212 s. The zero phase is arbitrary. 



light curve for the non-flaring time interval is plotted in Fig.|7] Con- 
trary to that obtained in 20 11, it shows a single broad peak. 

We finally performed a timing analysis of a series of 
RXTEfPCA archival observations carried out between 2009 April 
5 and April 15, obtaining a best fit spin period of 1209.4±1 s. The 
comparison of the pulse periods measured with XMM-Newton 
in 2004 and 201 1 and with RXTE in 2009 does not show signifi- 
cant spin-up or spin-down. 



4 DISCUSSION 

INTEGRAL showed that IGR J16418-4532 is an X-ray tran- 
sient with frequent flaring activity: it detected source outbursts 



above 10"^" erg cm"^ s"^ (20-100 keV) for about 1% of 
the INTEGRAL total exposure time of the source region 
l lDucci et al.l l2010h . This led to suggest IGR J16418-4532 as a 
member of the class of the supergiant fast X-ray transients (SFXTs, 
Sguera et al. (20q3); see|sidoli (2010) for an updated review of the 
SFXTs properties). It has been classified as a "candidate" SFXT be- 
cause the supergiant nature of its optical counterpart has not been 
confirmed yet, and the X-ray dynamic range was not as extreme 
as the prototypical SFXTs (typically, from 3 to 5 orders of mag- 
nitudes in X-ray intensity). Among SFXTs, IGR 116418^532 is, 
after IGR J16479~4514 and IGR J1848 3-0311, the source with 
the m ost frequent bright flaring activity jPucci et al.ll201(]l ISidolil 
I20T0I) . 

Our XMM -Newton observation allowed us to continuously 
follow the source variability with unprecedented detail and to re- 
veal, for the first time in this source, a wide dynamic range, typical 
of ot her so-called "inter mediate" SFXTs, like e.g. IGR J18483 - 
031 1 dSguera et al.ll2007h or IGR J16465-4507 jClark et alj2oTol) . 

The X-ray light curve observed in 20 11 showed two episodes 
of bright flaring activity (exceeding 10 counts s^^ in the first part 
of the exposure), separated by an interval of low intensity emission, 
where the source faded to less than ~0.1 counts s~^ (Fig. [8] mid- 
dle panel). The two flaring intervals (D and F) were characterized 
by different variability patterns: in the first one the flux reached the 
peak during the longest flare (~4000 s) and showed a variety of 
flare shapes, while in interval F the flares had a more regular ap- 
pearance, with a hint of a quasi-periodic pattern of a few ks, not 
related with the long spin period (see the lower panel of Fig. [8]l. 
A quasi-periodic flaring was also present in the central part of the 
2004 observation (time interval B in Fig. 1). As visible in the up- 
per panel of Fig. [8] the variability consisted of two sub-sets of four 
flares each, recurring every ~380-400 s. 

Quasi-periodic flares have been seen also in oth er HMXBs, 
includ ing some SFXTs (e.g. XTE J1739-302, see iDucci et al] 
I2OIOI) . For example, in the Be/X-ray transient EXO 2030+375 
flares recurring on timescales of 3.96±0.04 hr, and intensity oscil- 
lations with periods of 900-1200 s, were observed and explained 
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Figure 8. Close-up views of IGR J16418^532 light curves shown in Fig. 1. 
In the middle panel we show the 201 1 light curve (same as in Fig. 1) in log- 
arithmic scale to more cleai'ly show the high dynamic range now observed 
in IGR J16418^532. In the upper and lower panels two close-up views of 
the 2004 and 2011 observations are displayed, respectively, to better show 
the time intervals where quasi-periodic flaring activity is present. 



with the formation and disruption of a temporary and tran sient 
accretion disk around the neutron star jParmar et al.l 1 198^ . Al- 
though in EXO 2030-1-375 both the flare shapes (well described by 
a fast rise and an exponential decay) and the time recurrence were 
more regular, it is possible that a similar mechanism is responsi- 
ble fb£ttieJGRJ1641 8-4532 flaring activity. The model proposed 
by iTaam et al. predicts that similar recurrence time scales 

for the flares can be well reproduced by a relative velocity of the 
neutron star with resp ect to the outflow ing wind of about 1200- 
1400 km s~^ (see also lDucci et al.ll201(]| and references therein). 

We note a striking similarity between the IGR 116418^532 
light curve observed in 2011 and the prediction of hydrodynamic 
simulations for systems accreting in the so-called "transitional 
case" between "pure" wind a ccretion and "full" Roche lobe over- 
flow telondin & o"wei]| 19971) . According to these simulations, in 
this particular transitional accretion regime, when the mass donor 
is close to fill its Roche lobe, the HMXB can display an X-ray 
light curve very similar to what we observed in IGR J16418-4532, 
both for what concerns the variability time scale and its dynamic 
range (see Fig.jS] middle panel). The mass loss from the supergiant 
is dominated by the strong wind, but with the additional contribu- 
tion of a weak tidal gas stream, focussed towards the neutron star 
(hereafter, NS). This mechanism produces extreme variations in the 
mass accretion rate, mainly due to the dynamical interaction of the 
weak tidal gas stream w ith the accretion bow shock around the NS 
telondin & Owenlll997h . We are aware that these simulations have 
been performed in two dimensions (computing the gas flow within 
the orbital plane) and the real 3D case is much more complex, but 



nevertheless we think that the similarity with the observed light 
curve deserves attention and further investigation . 

Similarly, Bhattacharya & van den Heuvell ( Il99lh discussed 
the case of HMXBs with short orbital periods (P^3-4.5 days), 
such as IGR 116418^532, as sources where the orbits are typi- 
cally narrow making accretion via "beginning atmospheric Roche- 
lobe overflow" an attainable mechanism to produce X-ray emis- 
sion. In general, under the assumption that a supergiant star has a 
sharply defined radius (e.g. the photospheric radius), Roche lobe 
overflow starts as soon as such radius extends beyond the Roche 
lobe. However, in reality the supergiant star does not have a sharp 
edge or radius and above its photosphere there is still atmospheric 
material in the form of strong stellar wind. Consequently, already 
before the photospheric radius reaches the Roche lobe, a small 
part of the stellar wind can begin to flow towards the neutron star 
through the inner Lagrangian point. In principle, this could also 
favour the formation of a transient and temporary accretion disk 
whose disruption could produce the observed quasi-periodic flares 
from IGRJ16418-4532. 

Given the short orbital period of 3.7389 days, we plot in Fig.|9] 
the Roche lobe radius of the companion star at periastron, as a func- 
tion of its mass and for different system eccentricities (Eggleto^ 
Il983h . The radii of massive early type stars of different luminos- 
ity classes jVaccaetal.ll 19961) are also plotted for comparison. The 
32.8 kK best fit effe ctive temperatur^ of the candidate counterpart 
teahoui et al 12008^ indicates a 08.5-09.5 type star, which can be 
in a narrow circular orbit with the NS, avoiding Roche lobe over- 
flow. For example, a 09 la (08. 5 la) star with a m ass of 46 Mq 
(50 A/0) and a radius of 24 Rq jVacca et al.ll996l) in a circular or- 
bit is consistent with the optical and X-ray observations. Another 
viable possibility is a 09.5 III companion with a mass of 24 Mq 
and a radius of ~15 R©. In this case the eccentricity can range 
from e=0 to e'--^0.2 without exceeding the Roche lobe surface. We 
are aware that the above speculations on the optical counterpart are 
rather uncertain in the lack of an optical spectrum, but we here only 
note that an early type supergiant (and a SFXT nature) is compat- 
ible with the short orbital period. Moreover, a narrow orbit in a 
massive binary is required in the "transitional case" we suggested 
above to explain IGR 116418^532 properties. 

Comparison of the 2011 spectral results with the previous 
XMM-Newton observation shows a lower absorbing column 
density in 2011 than in 2004. However, in both occasions, the ab- 
sorption was well in excess with respect to the total Galactic col- 
umn density towards the source. It is possible that the variability of 
the absorption in the two observations is due to the different orbital 
phases covered by the two observations, but the uncertainties on the 
orbital parameters do not allow us to determine the absolute orbital 
phase covered by two spectra. The local absorbing column density 
due to the supergiant wind in a HMXBs is expected to vary along 
the orbit, even in case of a spherically symmetric wind, also in de- 
pendence of the inclination of the system. Further variability can 
be due to the presence of different gas structures (accretion wake 
around the NS together with tidal effects which induce a gas stream 
from the supergiant) which form becau se of the presence of the NS 
gravitational field (Blond in et aljl99"ll) . 

Another important result is that the possible hint for a soft 
excess present in 2004, is now well established by the 201 1 obser- 
vation and it is compatible with the presence of a ionized absorber. 
In particular, the cumulative X-ray spectrum during the flares ob- 



^ with very large uncertainties: Tdi can range from 36 kK to 10.6 kK 
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Figure 9. Roche lobe radius at periastron versus the companion mass, as- 
suming an orbital period of 3.7389 days, and differ ent eccentr icities. We 
oveiplot also the radii of early type stars, taken from Vacca et al.l fl996): di- 
amonds and triangles mark OB star with luminosity class III, their spectro- 
scopic and evolutionary masses, respectively (Table 6 in Vacca et al. 1996), 
vifhile asterisks and open squares mark supergiant stars (data are from Ta- 
ble 7 in Vacca et al. 1996), their spectroscopic and evolutionary masses, 
respectively, for different spectral types. 



The kind of X-ray variability (its higii dynamic range on tiie 
observed time scale), together with the candidate optical counter- 
part, the short orbital period and the quasi-periodic flaring activ- 
ity, all suggest that the X-ray emission from IGR J16418-4532 
is driven by a transitional accretion regime, intermediate between 
"pure" wind accretion and "full" Roche lobe overflow. 

In conclusion, we suggest here for the first time that this hy- 
pothesis could explain both IGR J16418-4532 X-ray behaviour 
and possibly other SFXTs with similarly short orbital periods. 
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served in 2011 (F) is well described introducing an additional 
ionized absorber, resulting in an absorber ionization state, 5, of 
125 I45 erg cm s~^, and an absorbing column density, Nnabsori, of 
15.8ti g X 10^^ cm"^. In the likely hypothesis of the wind photo- 
ionization, the above value of the ionization parameter ^ indicates, 
for example, the presenc e of highly ionized o xygen (O VIII) and 
neon (Ne X ions) (Kall man & McCravlll982h . whose lines could 
have been observed in the RGS energy range, if the source had 
been brighter and less absorbed. The average X-ray flux during 
flares, corrected for both interstellar and local (ionized) absorption 
is 5xl0~^^ erg cm^'^ s~^ (1-10 keV), which implies a ionizing 
X-ray luminosity L=10^'' erg s^^. Although the ^ value we ob- 
tained is very likely an average value, it is possible to calculate 
the distance, R, from the X-ray source of the main component of 
the absorbing ionized material. Using the measured Nnabsori value 
and assuming Nnabsori = nR, we obtain a distance R=5 x 10^" cm, 
which is compatible with the NS accretion radius, and well within 
the orbital separation (~10^^ cm) of the binary system. 



5 CONCLUSION 

The new XMM-Newton observations we have reported here 
allowed us to perform an in-depth investigation of the transient 
source IGR J164 18-4532. A high dynamic range has been observed 
for the first time, of about two orders of magnitude, leading to a 
firmer classification of IGR J16418-4532 as a member of the class 
of the SFXTs. 

We obtained a more precise value of the pulse period and we 
clearly established the presence of a soft X-ray excess in the flares 
X-ray spectrum, which we interpreted as due to the presence of 
ionized wind material, mainly located at a distance compatible with 
the NS accretion radius. The absorbing column density variations 
between 2004 and 201 1 are likely due to the different orbital phases 
covered between the two observations. The low intensity emission 
is both softer and more absorbed than during bright flares. 
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